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Abstract 

A Lorentz-covariant regularization scheme for effective field theories 
, with an arbitrary number of propagating heavy and light particles is 

given. This regularization scheme leaves the low-energy analytic structure 
£N| ' of Greens functions intact and preserves all the symmetries of the under- 

. lying Lagrangian. The power divergences of regularized loop integrals are 

controlled by the low-energy kinematic variables. Simple diagrammatic 
rules are derived for the regularization of arbitrary one-loop graphs and 
1 the generalization to higher loops is discussed. 
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■ 1 Introduction 



Effective Field Theory (EFT) is a consistent framework for calculating Greens 
functions below a certain scale Aeft when the fundamental theory is unsolvable. 
It is based on the principle that a physical amplitude can be calculated from the 
most general Lagrangian formulated with the relevant degrees of freedom and 
consistent with the symmetries of the fundamental theory. EFT then provides 
a systematic method of organizing an expansion in p/Aeft where p is a generic 
kinematic variable much smaller than Aeft- This expansion requires a regu- 
larization scheme which preserves the power counting in the small kinematic 
variables. 

EFT has been applied successfully to the calculation of Greens functions that 
involve strongly interacting particles at low-energies yj . One example is chiral 
perturbation theory (%PT) |^] where only quasi- Goldstone bosons of generic 

*dlehmann@jlab . org 
tprezeau@j lab . org 

■^Mail address: Thomas Jefferson National Accelerator Facility, 12000 Jefferson Avenue, 
M.S. 12H2, Newport News, VA 23606, USA 



1 



mass m are present and all the kinematic variables are small relative to A x sw 
lGeV. Dimensional regularization can be thought of as a natural regularization 
scheme for xPT, because dimensionally regularized loop integrals can only yield 
terms that involve powers of the small kinematic variables such as m/A x .Q 
Hence, dimensional regularization is a preferred regularization scheme because 
it allows the assignment of a chiral power to any Feynman graph through a 
naive dimensional analysis while preserving the symmetries of the problem. 

In the presence of heavy particles with generic masses denoted by M > Aeft, 
naive dimensional analysis is no longer valid because the loop integrals also yield 
terms of the form M/Aeft which spoil power counting S. 

For a single heavy particle at low energies, the dependence on M can be 
removed from the propagators by expanding the EFT Lagrangian in inverse 
powers of the heavy mass as is done in heavy quark effective field theory and 
heavy baryon chiral perturbation theory jl. The 1/M expansion yields loop 
integrals that satisfy low-energy power counting. However, the 1/M expansion 
does not always reproduce the correct low-energy analytic structure. An exam- 
ple of this is the scalar form- factor of the nucleon where the 1/M expansion fails 
to reproduce the anomalous threshold 0. Therefore, the 1/M expansion near 
the anomalous threshold will fail to converge and an infinite number of terms 
must be re-summed. This re-summation restores manifest Lorentz-invariance. 
The retention of manifest Lorentz-invariance from the beginning avoids these 
problems. 

A manifestly Lorentz-invariant regularization scheme which preserves the 
low-energy analytic structure without spoiling power counting is desirable. In 
the presence of a single heavy particle, such a regularization scheme called in- 
frared regularization was devised for one-loop Feynman graphs by Becher and 
Leutwyler || building on work by Ellis and Tang |Q. The underlying idea 
of infrared regularization is the separation of power-counting violating terms 
from the dimensionally regularized loop integrals, and their absorption into the 
low-energy constants (LECs) of the EFT. The power-counting violating terms 
are proportional to fractional powers of M which involve non-integer values of 
the dimension d. For example, for the nucleon self-energy graph of Fig. |l|, the 
infrared regularization separates out terms of the form M d ~ A /(d — 4). This 
separation is achieved by extending the integration domain of the Feynman pa- 
rameters appearing in Greens functions. Inspired by this work, a regularization 
of one-loop graphs with two propagating heavy particles was recently developed 
in Ref. || and referred to there as EFT dimensional regularization (EFTDR). 

In this paper, a consistent, Lorentz-invariant, natural regularization scheme 
is presented which generalizes EFTDR to an arbitrary number of heavy and 
light particles and to higher loops. This regularization systematically separates 
out all terms involving fractional powers of M implementing a low-energy power 
counting scheme. To one-loop order, it is proven that the low-energy analytic 
structure of the Greens functions is left intact by this procedure such that the 
regularization is valid throughout the low-energy region. It is demonstrated in 

1 In contrast to cutoff regularization which introduces powers of the cutoff. 
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Figure 1: Contribution to the nucleon self energy. 



the appendix that this general formalism reproduces the one and two heavy 
particle sectors as first derived in Refs. |6[ |). 

In Sect. the analytic properties of Feynman graphs are reviewed employing 
the Landau equations. It is shown that the singularities of Feynman integrals 
can be classified according to subgraphs of the original graph. In Sect. ||, the 
regularization scheme is introduced. It is shown that the fractional powers of M 
can be associated with subgraphs of the original graph and can be systematically 
separated out. To one- loop, simple diagrammatic rules are given to achieve this 
separation. The application of the formalism to multi-loop graphs is discussed. 
The results are summarized in Sect. 0. 



2 Analytic structure of Feynman graphs 

Consider a Feynman graph G with / scalar internal lines and L loops in d 
space-time dimensions, 

where /i is the renormalization scale of dimensional regularization. The momen- 
tum of each internal line i is a linear combination of the loop momenta 
and the external momenta Ptf . The Feynman graph thus depends on the parti- 
cle masses wii and a minimal set of scalar combinations P v ■ P v > of the external 
momenta. In a general EFT, fermion propagators and derivative couplings give 
rise to a momentum-dependent numerator in Eq. (|l|). Since the analytic prop- 
erties of a Feynman graph only depend on the structure of the denominator as 
will be shown below, momentum dependent numerators can be neglected for the 
purpose of this paper. This simplification is not a limitation of the formalism 
as Feynman graphs involving momentum dependent numerators can always be 
decomposed into a combination of integrals of the form given by Eq. (|l|). 
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The propagators of Eq. (|l]) can be combined with Feynman integrals to yield 

/ 

I . L Jd7 S ( 1 ~ E <*j) 

d a k e 3=1 



Iq = i J r(J) M i ( 4 - rf ) r f[ \da k 8(a k )] f f[ 
J -°°k=i J 1=1 



i=i y ' 



E -m?)+ie 

i=l 



(2) 



where the ^-functions have been introduced for later convenience. In general, 
singularities of Iq will arise when the contour of integration gets pinched be- 
tween two or more poles of the integrand (pinch singularities) or between poles 
of the integrand and endpoints of the integration contour (endpoint singulari- 
ties). In either case, the integration contour cannot be distorted away from the 
pole of the integrand resulting in a singularity of the integral itself. 

A useful tool for investigating analytic properties of Feynman graphs is given 
by the set of Landau equations.^ They constitute a necessary but not sufficient 
condition for the occurrence of singularities in Iq and their solution determines 
the locations of the singularities in the space of masses and external momenta 
as well as their locations in the space of integration variables a*. In the rep- 
resentation (J2j) of the Feynman graph, the Landau equations for Iq are given 

by 

for each i: either a, = or <jf = rrq , (3) 
for each I, fi: ^2^-^ = 0, (4) 

where the c^'s must be normalized to satisfy the 5-function constraint, 

X> = 1. (5) 

i=l 

An immediate consequence of Eq. (m is the vanishing of the denominator in 
Eq- @, E[=i a i(li ~ m f) = • Eq. means that an internal particle i of the 
Feynman graph G is either on-shell or the Feynman parameter on of line i is 
zero. If oti — 0, line i does not contribute to the remaining Landau equations 
(g). In this case, the set of Landau equations for G is identical to the Landau 
equations for the subgraph g of G that is obtained by contracting line i to a 
point ^] 

A solution to the Landau equations where all the o^-'s are non-zero, i.e., all 
internal lines are on-shell, is called the leading singularity of G; singularities with 
n internal lines off-shell are called ( sub ) n -leading singularities. It follows that 

2 See Refs. ^ |ic|] for a review. 

3 It is important to emphasize that the subgraphs thus obtained are not Feynman graphs 
and do not represent physical amplitudes. Henceforth, these subgraphs will be represented by 
lower-case letters. 
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every (sub)™-leading singularity of a Feynman graph G can be represented by 
the leading singularity of a subgraph g where the corresponding n internal lines 
of the original graph have been contracted. The complete set of singularities 
of a Feynman graph G is then given by the leading singularities of each of its 
subgraphs. 

For one-loop graphs, a more practical formulation of the Landau equations 
can be derived. Routing all internal momenta to flow around the loop, Eq. (||) 
takes on the form 

for each i with on ^ 0: ct^ = , (6) 

i=i 

where only uncontracted lines j contribute. Multiplying this equation by 
for each uncontracted line i yields a set of linear equations 

■ 1j) <*j = ■ (7) 
Using the on-shell condition (^), the matrix (g, ■ qj) can be re- written as 

Oij = ■Qj) = — . ( 8 ) 

which is now independent of the loop momentum The location of the sin- 
gularity in the space of masses and external momenta is determined from the 
condition 

det % = , (9) 

while the corresponding location in Feynman parameter space is obtained by 
solving the zero-eigenvalue problem ([?]) and normalizing the ctj's according to 
Eq. (||) . A singularity does not arise on the physical sheet when the solution in 
parameter space lies outside the domain of integration of the Feynman graph. 



2.1 An example: the triangle graph 

As an illustration consider the triangle graph shown in Fig. ||(a) where the solid 
lines represent heavy particles of mass M and the dashed line represents a light 
particle of mass m. The corresponding subgraphs are shown in Figs, ^(b)-(g). 

For the triangle graph, the 5-function appearing in Eq. (||) constrains the 
domain of integration to a plane in the (a±, a%, a3)-space. Fig. |^ shows the 
integration domain and the locations of the singularities in this plane. The 
integration domain (light shaded area) is bounded by the lines on — 0. A (sub) 1 - 
lcading singularity corresponding to the subgraph where line i is contracted lies 
on the line on = 0. A (sub) 2 -leading singularity represented by the subgraphs 
where lines i and j are contracted lies at the intersection of the lines cti — 
and ctj = which is a corner of the integration domain. Finally, the Feynman 
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(e) (f) (g) 

Figure 2: Triangle graph (a) and all its subgraphs (b)-(g). 



parameters that correspond to the leading singularity form a surface given by 
the dark shaded area in Fig. [| which penetrates the integration boundaries at the 
location of the (sub) ^leading singularities. In Feynman parameter space, the 
leading singularity is a pinch singularity while all the sub-leading singularities 
occurring on the physical sheet lie on the integration boundary and are therefore 
endpoint singularities. 

It is instructive to solve the Landau equations for each subgraph: 

• (Sub) 2 -leading singularities. For the heavy-particle tadpole subgraph 
Fig. g(e) with lines 2 and 3 contracted, Eqs. (fjj) and (g) yield immedi- 
ately 

eti ^ , a 2 = a 3 = , 
ai M 2 = => M 2 = . 

Similarly, the light-particle tadpole Fig. ||(g) yields a singularity at m 2 = 0. 
In Fig. |[ the heavy-particle tadpole singularities are denoted by ®, the 
light-particle tadpole by ■. 

• (Sub) 1 -leading singularities. For the self-energy subgraph in Fig. §(b) 



G 



Figure 3: Integration domain and location of singularities of the triangle graph. 



condition (|^) reads 
ai = , a 2 ^ , a 3 ^ , 

,1 M 2 +rn 2 -P 2 



M 2 

M 2 +rn 2 -P? 



m 



= [Pi - (M + m) 2 ] [P 2 2 - (M - m) 2 ] = . 

(11) 



The corresponding zero-eigenvectors are easily determined. The singular- 
ities are located at 



(M - 



Pi = (M - mf 



for a 2 



(7? 



M 



for 



M + m 
m 



, a 3 



a 2 = 



M + m 
M 

, "3 = 



(12) 
(13) 



M — m M — m 

Solution (12) (denoted by • in Fig. ||) is recognized as the physical heavy- 
light-threshold, while the pseudo-threshold ( J13| ) (denoted by o) does not 
appear as a singularity on the physical sheet since it lies outside {a 2 < 
0) the domain of integration of the Fcynman graph (|J). Similarly, the 
thresholds for the bubble subgraph in Fig. ||(d) are 



P 2 = AM 2 



P 2 = 



for 



a i = ot 2 = 



a 3 = , 



not on the physical sheet 



(x in Fig. |); (14) 
(15) 
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• Leading Singularity. The leading singularity of the triangle graph, corre- 
sponding to oti 7^ (i = 1, 2, 3), is a manifold rather than a point in the 



space of independent external momenta P 2 , Pf , and P 2 
satisfies 



M 2 

2M 2 -P 2 
M 2 +m 2 -l 



2AI- 



M 2 

M 2 +m 2 -P? 



M 



(P 1 +P 2 ) 2 and 



(16) 



To this manifold in momentum space corresponds a manifold in Feynman 
parameter space given by the dark shaded area in Fig. 0. 



3 EFT Dimensional Regularization 

In this section, a natural regularization scheme (EFTDR) for an EFT applicable 
below a scale Aeft is established. It is the low-energy kinematic variables that 
control the divergences of the loop integrals regularized using EFTDR. The 
low-energy domain is defined by the relations 

|P„-P/-M 2 |«A| FT , 
P v ■ Pv' , Pv ■ Pv' , m 2 < A| FT , (17) 
M > Aeft , 

where P£ and p% are the heavy and light particle external momenta respectively, 
while M and m are generic heavy and light particle masses respectively. 

EFTDR is implemented by consistently separating a loop integral Ig corre- 
sponding to a graph G into two parts, 

I G = 1 G + R G , (18) 

where the entire low-energy analytic structure of Iq is collected into Iq while 
all the terms multiplied by factors of the form M ed ~ n (referred to as fractional 
powers of M for fractional space-time dimension d) where I and n are integers, 
are collected into Rq, the regular partJ^ The analytic structure of Iq and Iq are 
identical in the low-energy domain such that Rg can be accounted for through 
a renormalization of the low-energy constants of the EFT. The separation ( p"8| ) 
will be achieved by systematically separating out the subgraphs representing 
singularities that lie outside of the low-energy region; these subgraphs will be 
referred to as regular subgraphs. 

In the example of the triangle graph, the only singularities that lie outside of 
the low-energy region are the M — singularities represented by Figs. |](e)-(f). 

4 Tho normalization condition (^|) cannot be fulfilled. 

5 The terminology is taken from Ref. |f and refers to the fact that to one-loop, Rq is 
analytic in the low-energy region. 
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Figure 4: A regular subgraph composed of heavy particle propagators and low- 
energy momentum insertions represented by wiggly lines. 



Since a factor M will be singular at M = for appropriate d, separating 
out the M — - singularity will ensure that the terms multiplied by fractional 
powers of M are separated out. 

Generally, subgraphs containing only heavy particle propagators and low- 
energy momentum insertions are regular subgraphs that give rise to fractional 
powers of M. Indeed, the Landau equations (^) for the regular subgraph shown 
in Fig. |] are 



M 



1 i 



7=1 



M = . 



(19) 



Therefore, these subgraphs cannot have singularities in the low-energy region 
and belong in Rq. 

The subgraphs that lead to fractional powers of M can be isolated by ex- 
tending the Feynman parameter domain of integration in Eq. where the 
0(ak) define the integral boundaries. The product of 0's can be rewritten: 



*(i-E>*)II[i 



k=l 



i=l 



1 r 1 

fc=i ^ t=i 



E [-ei-oHE-Bi-oj)] 



»)j=i 



(20) 



It is noted that each 6{—a.i) contains the point on = which, as discussed in 
the previous section, represents a contracted line in G that leads to a subgraph 
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of G. In light of this, each term in Eq. ( |20| ) defines a domain of integration A 9 , 



\ = S(l-^2a k )l[[-9(-a i )] 



(21) 



k=l 



which can be associated with a particular subgraph g of Gjj Rq is then defined 

as 



r g = i I r(i) f i L ^- d '> 

where 



k=l 



1 d d h r ' 



j (2tt)° 



1 mj)+ie 



= E A.. 

regular g's 



(22) 



(23) 



It is instructive to discuss the triangle graph / 3 as an illustration of the method. 
Tab. [I] shows the association of each term in the ^-expansion with the corre- 
sponding subgraph and domain of integration. As mentioned above, the only 
regular subgraphs are the heavy particle tadpoles located at (0.1,0.2,01.3) = 
(1, 0, 0) and (0, 1, 0) in Fig [| The integration domains in Eq. ( |20| ) that contain 
these points as endpoints are: 8(— a2)8(— a 3 ) and 9(— a\)8(— a 3 ). These inte- 
gration domains do not include any of the low-energy singularities; this implies 
that i?3 is analytic in the low-energy domain. -Z3 is then given by 



la — I?, — R?, 



(24) 



and possesses the same low-energy analytic structure in the small kinematical 
variables as J3 but does not contain any fractional powers of M. This is evident 
by looking at the integration domains corresponding to Eq. ( p4[ ) and drawn in 
Fig. I 









\ ) 
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Figure 5: Decomposition of the integration contour for the triangle graph 



For a general Feynman graph G, a particular domain of integration A g 
associated with a subgraph g will never select a singularity associated with any 
subgraph of g. Indeed, the domain of integration A g will never have as an 

6 The subscript g represents a set of indices labeling the internal lines of the subgraph g 
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subgraph g 



domain function A 



integration domain 



bJO 

CO 




Pi 



s(i - E «i) 





cd 
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• ► P 2 



-^(-01) 

x<5(l- E«, 

i=l 



Pi 




\ 



-0(-a 2 ) 



Pi 



x<5(l - E a t) 

1=1 



P 2 




P 



-0(-a 3 ) 






-ai) 0(-a 2 ) 

3 

X(5(l - E a ») 

1=1 




be 




(-ai) 0(-a 3 ) 

3 

X(5(l - E ai) 

8 = 1 




0(-a 2 )0(-a 3 ) 

3 

X(5(l - E ai) 

i=l 




1 


<*- 







Table 1: Integration domains for the triangle graph 
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endpoint the location in a-space of any singularity associated with a subgraph 
of g. It follows that a subgraph g that contributes to Iq can never give rise to 
fractional powers of M since they can only arise from regular subgraphs that 
are subgraphs of ~g. This proves that Rg contains all the fractional powers of 
M. To one loop, it is also straightforward to show that Rg is analytic in the 
low-energy region. 



3.1 One-loop graphs 

For one-loop graphs, more explicit results can be obtained. Evidently, any 
subgraph of a regular subgraph, obtained by contracting out further lines, is 
also regular. This gives rise to the notion of a minimally contracted regular 
(MCR) subgraph, i.e., a regular subgraph g that is not a subgraph of any other 
regular subgraph with more lines than g. In EFTs to one-loop order, there 
are at most two different MCR subgraphs of a given graph and every internal 
heavy-particle line belongs to exactly one of them (cf. App. |a|). As an example, 
the two tadpoles of Figs. 0(e)-(f) are the MCR subgraphs of the triangle graph. 

Since by construction different MCR subgraphs do not have common lines, 
they cannot have subgraphs in common. Any regular subgraph g is thus a 
subgraph of exactly one MCR subgraph g. The domain functions (HJ) of a 
MCR subgraph g and all its subgraphs g C g can be re-summed to yield 

E A g = 5(1 - E «*) {l - E 9 ^ + ' ' ' } II ( 25 ) 

(26) 



gQs 



k=l 
I 



'to 



fc=l ie$ j^g 

The domain function T>r g introduced in Eq. (|2^) is then 
and the full regular part is given as 

t>oo />0 £ 

R G = i I T(I)^- d ^2 II da ' / IK-^^-E"*) 

a Jo i£B j&B k=l 

x E a <( 9 « -TO <) +ie 

^ ' L z— 1 

The momentum integration can be carried out explicitly |l0| yielding 

,o I 



# G = ^ e / n da * / n(- d ^) ^ - e t c - k ] 



. id/2-I 



(27) 



(28) 



(29) 
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where: 



C = a * n 



+ m 2 - ( qi - qjf 



(30) 



The matrix flij was defined in Eq. (0) and the factors have been collected in 



«/ = (-) 



r(/-d/2) 



(31) 



16tt 2 (47T^ 2 ) d / 2 - 2 ' 

A useful representation of Rg with a single non-compact integral can be 
derived by inserting the identity 



/oo 
dA5(A-^a fc ) 
fcgfl 

into Eq. ( |29| ) and performing the change of variables 

a, = (1 - A)zi i ^ g . 

This yields for i?^: 

/oo 
dAAlol-Hl-A)'-! !- 1 
» - 

x /' ( n ) *(! - e **) /' ( n ^ ) ^ - e **) 

(1 — A) 2 S ' ZiftijZj + A 2 ^ ] ZiilijZj 



(32) 



(33) 



2A(1 — A) ZiflijZj — ie 



d/2-I 



(34) 



where |g| denotes the number of internal lines belonging to the MCR subgraph 
0- 

In this representation, it is straightforward to prove that Rq is analytic in the 
low-energy region. As stated above, every internal heavy-particle line belongs 
to one of the two MCR subgraphs 91 or Q2^\ Denoting the set of internal light- 
particle lines by I, the elements of the matrix appearing in Eq. (|3^) are (cf. 
Eq. © in App. 

(0( P 2 ) ijel 

}0(p)M j$l 



M 2 + C(p 2 ) y'ejh or i,jeg 2 

{-M 2 + 0(p)M ieji, jeg 2 



7 Here g2 is understood to be empty in cases where there is only one MCR subgraph. 
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yielding for the denominator of Eq. (|34j) with Q — Qi- 
i 



aitoijotj = (1 - A) 2 Y ZiM 2 Zj - 2A(1 - A) ^ ^M 1 

A 2 £ z 4 M 2 z, + 0(p)M (36) 



M 2 (1 + ^) 2 A 



2 \2 



1 



(l + i/)A 



0(p)M, 



where z/ = ^ with < v < 1 and A > 1. Thus, the denominator never 

vanishes proving that Rq is always analytic for one-loop graphs with arbitrary 
numbers of heavy and light particle propagators. The integral over A can be 
performed by expanding the square bracket in Eq. (|36| ) . To leading order in this 
expansion, this gives 



(_l)/-|B|-l M d-2/ K/ 



c (d-/-i)r(| fl |) 



vf- 11 + MO{ P ) , (37) 



where the A-integral was performed at d sufficiently small to drop the upper 
boundary; d can now be analytically continued to 4 space-time dimensions. 
Eq. (|37| ) can be used to calculate to leading order the regular part of R3 keeping 
in mind that there are two MCR subgraphs for the triangle graph: 

_ 1 r(3-d/2) 
3 16^ 2 (4^ 2 ) d / 2 - 2 (d-4)(d-5) 1 j ' K > 

Eq. ( |38| ) agrees with the result given in Ref. || . 

Eq. ( |37| ) explicitly shows that Rq is proportional to a fractional power of M. 
This implies that Rg and Iq transform separately under the symmetry trans- 
formations that leave the Lagrangian invariant [||. This ensures that EFTDR 
preserves all the symmetries of the underlying Lagrangian. For general one-loop 
graphs, Eq. ( |37| ) also shows that the non-compact integrals appearing in Iq can 
be performed quite generally. 

The Taylor series expansion of Rg in powers of p is truncated at the order 
of the EFT Lagrangian. It follows that the Taylor-polynomial of Rq can be 
absorbed into the low-energy constants of the EFT Lagrangian. Once Rq has 
been used to define the domain of integration of Ig, it need not be considered 
further. 

For a general one-loop Feynman graph G, the EFT dimensional regulariza- 
tion scheme is summarized by the following simple rules: 
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1. Label all of the internal lines i of the graph G by the corresponding Feyn- 
man parameter on. 

2. Determine the minimally contracted regular subgraph(s) g of G. 

3. For each minimally contracted regular subgraph g, integrate the Feynman 
parameters corresponding to internal lines of g from to +oo; integrate 
the remaining Feynman parameters from — oo to inserting a factor of -1 
for every integral. 

4. Sum these contributions for all minimally contracted subgraphs to obtain 
the regular part Rq. 

5. The EFT dimensionally regularized Feynman integral is given as 

Ig = Ig — Rg ■ 

3.2 Discussion of multi-loop graphs 

It is emphasized that the regularization scheme developed above is not lim- 
ited to one-loop graphs. Indeed, the Landau equations are valid to arbitrary 
number of loops and allow a determination of those subgraphs that give rise to 
fractional powers of M. The contributions from these and only these subgraphs 
are collected in Ro and need to be separated from the dimensionally regularized 
Feynman integral Iq to obtain the EFT dimensionally regularized integral Iq. 
The separation 

I G = 1 G + Rg (39) 

is again achieved by extending the domain of integration as explained above. 
As a novel feature at two or more loops, subgraphs that give rise to fractional 
powers of M may now contain light internal lines and may give rise to low- 
energy singularities. An example of such a subgraph is a tadpole attached to a 
triangle graph as shown in Fig. ^|(a): the tadpole part of the graph gives rise to 
fractional powers of M, while the triangle part contains low-energy singularities. 
The appearance of these graphs does not constitute a problem since their low- 
energy singularities are not new types of singularities but rather originate from a 
lower number of loops. Even after separating out Rg, all the thresholds and the 
associated analytic structure of these low-energy singularities are still contained 
in the low-energy part Iq through subgraphs with fewer loops. In the two-loop 
example shown in Fig. |6| the low-energy singularities are identical to those of the 
one-loop triangle graph which are retained in Iq according to the regularization 
scheme and the contribution of Fig. ^| corresponds to the triangle graph with a 
renormalized four- heavy-particle vertex shown in Fig. ||(b) . 
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(a) (b) 

Figure 6: Figure (a) shows an example of an unnested two-loop graph that gives 
rise to fractional powers of M while containing low-energy singularities. Fig- 
ure (b) shows the triangle graph with renormalized four-heavy-particle vertex. 



4 Summary 

In summary, a manifestly Lorentz-invariant regularization scheme which does 
not violate power counting was derived. This scheme separates a Feynman in- 
tegral Iq represented by a Feynman graph G into two parts, Iq — Ig + Rg, 
and is referred to as EFTDR; the entire low-energy analytic structure of Iq is 
collected into Iq while all the low-energy power counting violating contribu- 
tions are collected into the regular part Rg- This separation is achieved by 
identifying the power counting violating contributions to particular solutions of 
the Landau equations which are represented by subgraphs of G. In EFTDR, 
each subgraph g of G is associated with a domain of integration A g in the space 
of Feynman parameters such that the power counting violating terms are sep- 
arated out by extending the domain of integration. Simple diagrammatic rules 
were given which allow the construction of Iq (and Rg) for one-loop diagrams 
with an arbitrary number of heavy propagating particles. Rg can be accounted 
through a renormalization of the LECs of the EFT Lagrangian and discarded; 
the EFTDR regularized Feynman graph Iq is then Iq. It was shown that Rg 
and Iq transform separately under the transformations that leave the EFT La- 
grangian invariant such that all the Ward identities are satisfied by the EFTDR 
regularized Greens functions of the EFT. The application of EFTDR to higher 
loops was discussed and it was emphasized that the power counting violating 
contributions can still be separated out systematically by extending the domain 
of integration in the space of Feynman parameters. EFTDR was applied to 
the one-heavy particle and two-heavy particle sectors to one loop and shown to 
reproduce the results in Refs. g|. 
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Figure 7: Cutting a one-loop graph 
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A Minimally contracted regular subgraphs 

Consider a general one-loop graph G formed with heavy (mass M) and light 
(mass to) particles and external momenta P£. By cutting two different lines i 
and j the graph is divided into two connected parts and , see Fig. ^. 
Four-momentum conservation at each vertex of G implies that the squared mo- 
mentum flowing through the cut, (g; — qj) 2 , equals the squared total external 
momentum attached to either one of the two pieces, 

(li 13? = (£(±)P„) 2 G [«] = (-^(±)P,0 2 G [«] , (40) 

V v' 

where the momenta of internal lines are routed to flow around the loop 
and the sign of the external momenta Ptf depends on whether the momenta 
are incoming or outgoing. Since in an EFT all external three-momenta are 
0(p 2 ) <§C Aeft and creation/annihilation of heavy particles is excluded, it fol- 
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(a) (b) 
Figure 8: Box graph (a) and cross box graph (b) 



lows immediately that^j 



'0(p 2 ) i, j light 

M 2 + 0(p)M i light, j heavy 

(2M) 2 + 0(p)M i, j heavy 

P(p 2 ) i, j heavy 



(41) 



An example for the case (^T|a) is found when cutting the two heavy-particle lines 
of the box graph Fig. ||(a), while the crossed box graph Fig. ||(b) corresponds 
to case (fy]b). 

The crucial observation is now that for any three heavy-particle internal lines i, 
j, and k: 



(ft - q 3 ) 2 = 0(p 2 ) 
and 

(q, - q k ) 2 = 0(p 2 ) , 



{q, - q k f = 0(p 2 ) 



(42) 



This transitivity relation follows from Eq. 



and 



(- V(±)px l «, = (V(±)p„) G «> + (E(±)^) 



(43) 



see Fig. ^|. Any two heavy-particle lines i and j of a graph G are said to belong 
to the same minimally contracted regular (MCR) subgraph g if and only if they 
satisfy condition Eq. (|d]b), 



i, i G 



- <z,) 2 = C7(p 2 ) 



(44) 



This defines the MCR subgraphs of G completely. In light of Eq. (|42j), differ- 
ent MCR subgraphs cannot have any internal lines in common. On the other 
hand, every heavy-particle internal line belongs at least to one regular sub- 
graph, namely the tadpole graph corresponding to that line. Therefore, the set 
of heavy-particle internal lines of a graph G can be split into disjoint subsets; 
each subset is diagrammatically represented by a subgraph g of G that consists 

8 It is to be noted that Eq. (|4l[) remains true for all subgraphs of G. Indeed, the l.h.s. is 
unaffected by the contraction lines to a point. 
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Figure 9: Transitivity of cuts 



of all the lines belonging to the subset and of no other lines. These subgraphs 
are the MCR subgraphs. They do not have low-energy singularities, since by 
construction all (effective) external momenta are small while all internal lines 
are heavy. 

In the box graph Fig. ||(a) the two heavy-particle lines 1 and 2 do not belong 
to the same MCR subgraph; in this case the tadpole graphs formed with lines 
1 and 2 are both MCR subgraphs. In the cross box graph Fig. |8|(b), however, 
there is only one MCR subgraph, namely the bubble graph formed with lines 1 
and 2. 

It is noted that at most two different MCR subgraphs of a given one-loop 
graph exist. This can be seen as follows: assume the existence of more than 
two different MCR subgraphs and choose three heavy internal lines i, j, and k 
belonging to different MCR subgraphs. It is easily verified that 

(<ft -q 3 f = AM 2 + 0{p)M\ 

and I ( qi - qk f ^ AM 2 + <D(p)M . (45) 

- Qkf = 4Af 2 + 0(p)M J 

However, since (qi — qk) 2 must satisfy Eq. (f4l|), lines i and k have to belong to 
the same MCR subgraph, 

( qi - q k ) 2 = 0(p 2 ) , (46) 
in contradiction to the above assumption. 

B One-loop examples 

This formalism reproduces the prescriptions given for the one-loop graphs in 
Figs. |l^(a) and |ll](a) first derived in Refs. || 0]. Fig. |l^(a) is the most general 
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(a) (b) 

Figure 10: (a) General one-loop graph with one heavy particle, (b) The corre- 
sponding minimally contracted regular subgraph. 



one-loop graph with a single heavy particle line running through it; there are 
H — 1 soft momentum insertions entering the heavy particle internal line and 
I — 1 soft momentum insertions entering the light particle internal line. 

The regular part for the graph in Fig. |l(j(a) corresponds to the MCR sub- 
graph in Fig. ^(b) where all the light internal lines have been contracted. 
Eq. (H]) gives 

/oo 
dA A^-^l-A)'- 1 

which is the result first obtained by Becher and Leutwyler || for a one-loop 
graph with I light particle propagators labeled 1, . . . ,1 and H heavy particle 
propagators labeled I + 1, . . . , I where I — I + H. Indeed, Eq. (47) describes a 
Feynman parameterization where all the light particle propagators were com- 
bined into a single denominator, all the heavy particle propagators were com- 
bined into another denominator, and these two denominators were subsequently 
combined with the Feynman parameter A integrated between one and infinity; 
this is precisely the regular part in Bcchcr-Lcutwylcr's infrared regularization. 

For a Feynman graph with two heavy particle lines running through it, 
Figs. [ll](a)-(b) are the most general diagrams that can be drawn.[] For the 
Feynman graph shown in Fig. |ll|(a), there are two MCR subgraphs obtained 
by contracting all the light particle internal lines as well as one of the heavy 
internal line running through the graph. The regular part is given by the sum 
of the MCR subgraphs. Assuming I light particle internal lines and that the top 
heavy particle internal line is segmented by hi — 1 soft momentum insertions 

9 For example, the triangle graph is obtained from Fig. |ll|(a) by removing one of the light 
particle internal lines that runs through the graph. 
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(a) (b) 
Figure 11: (a) Generalized box graph, (b) Generalized cross box graph. 



while the bottom heavy particle internal line is segmented by /12 — 1 soft mo- 
mentum insertions, a change of variables similar to Eq. (|3^) can be performed 
on Eq. ([54]) to cast it in the form in which it appears in Rcf. || . For the MCR 
obtained by contracting the h n (n = 1, 2) propagators of one of the heavy lines, 
inserting the integrals 

d\df3 6(\- ot k )8{p-- Y, afc ) ( 48 ) 

k=l + l k=l+h„ + l 

and performing the change of variables 

a.i = (1 - X)z l i = l,...,l 
a i+ i = A(l - f$)xi i — l,...,h n (49) 



Ui+i+h„. = X/3yi i = 1, . . . , H - h r 



gives 



/>oo />oo 

R= V -ki / dA / df3X H - 1 {l^X) l - 1 f] H - h "- 1 (l-f3) 



hn-1 



n=l,2 




H-h„ 



J] dyA [l[dz k 

3 = 1 J \fc=l / 



H-h n h„ 



2/6)*(l -$>c)[C-ie] rf / 2 -', (50) 



a=l b=l 



where H = hi + hi. This is just the regular part given in M. The generalized 
cross box graph in Fig. |ll|(b) does not have a two-heavy particle threshold, 
and its regular part is constructed from a single MCR subgraph obtained by 
contracting all the light internal lines. 
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